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1. INTRODUCTION 

Surface acoustic wave (SAW) devices have numerous applications in the fields of actuators and 
sensors, including the microelectromechanical systems (MEMS) actuator [1]. SAW sensors have been 
applied in various fields, such as mechanical, biological, chemical, gas, and microfluidics, to measure 
physical quantities including temperature, pressure, viscosity, mass, stress, current, gas composition, 
magnetic field, and the spatial localization of partial discharges [2]-[9]. Recently, the wireless magnetic 
sensors with low energy consumption are essential for developments in the field of the internet of things, 
autonomous vehicles and the factory of the future. Common constructions for magnetic SAW sensors usually 
have a magnetic sensitive materials, an intermediate layers, an electrodes and a piezoelectric substrates 
(magnetic sensitive materials/intermediate layers/electrodes/piezoelectric substrates). Several recent studies 
indicate that common magnetostrictive sensitive materials, including Ni, FeGa, FeCo, and FeNi, have been 
combined with quartz and LiNbO; piezoelectric substrates to measure magnetic field intensity and current with 
high sensitivity [4]-[6], [10]-[17]. The results show that the thickness of the magnetic materials significantly 
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affect the sensitivity of the sensor [4], [11], [15]. Notely that the simulation approach was used to optimize 
the design parameters of sensors, in which the electromechanical equations set are solved numerically 
through finding boundary conditions to satisfy the CHRISTOFFEL equation and using perturbation theory 
[11]. In fact, it is difficult to determine the boundary conditions, and thus consequently this approach is not 
convenient to implement. Currently, the most popular simulation method for SAW sensors is using finite 
element method (FEM) through commercial software tools such as COMSOL and ANSYS. The survey 
shows that many studies have used FEM to simulate gas sensors, pressure sensors, magnetic field sensors, 
and SAW filters [6], [18]-[23]. Namely, the FEM with the COMSOL tool was used to evaluate the influence 
of the intermediate layer on the sensor properties of the magnetic SAW sensor [16]. 

Thus, in this work, we used the FEM with ANSYS software to simulate the SAW-MO sensor with 
FeNi/IDT/AIN structures. In this SAW-MO sensor, FeNi is selected because it has good sensitivity, low 
nonlinearity and small hysteresis; AIN is a piezoelectric material and has a large surface wave velocity that 
leads to being able to achieve the high sensitivity. The thickness of the FeNi and AIN layers changed during 
the simulation to determine the optimal parameters for the sensor structure to achieve the highest sensitivity. 

The rest of the paper is organized as follows: section 2 presents the simulation process of the paper, 
which is divided into smaller sections, whereby section 2.1 explains the principle and structure parameters of 
the sensor, while section 2.2 gives mathematical model for FEM simulation, and section 2.3 presents FEM 
simulation process. Section 3 shows the results and discussions of this research paper and, finally, section 4 
displays the conclusive findings for this research paper. 


2. SIMULATION PROCESS 
2.1. Principle and structure parameters of the sensor 

Figure 1 illustrates the SAW-MO sensor structure. The sensor comprises an AIN piezoelectric 
substrate, which has two sets of electrodes, namely the input (IDT-in) and output (IDT-out) electrodes. These 
electrodes sandwich the FeNi sensitive layer (in a delay-line form). When a voltage is applied to the IDT-in, 
the piezoelectric substrate generates a surface acoustic wave through the reverse piezoelectric effect, which 
propagates through the FeNi sensitive layer to the IDT-out, where a forward piezoelectric effect occurs and 
forms a voltage on the IDT-out. The interaction between the measured magnetic field intensity (H) and the 
FeNi magnetic sensitive layer changes the frequency of this signal (or the SAW velocity). The structural 
parameters of the sensor are detailed in Figure 2 and Table 1. In the simulations, the sensor has a constant 
length (in the x, axis) and depth (in the x3 axis), and initially, the sensor is simulated with d = 10 um 
(corresponding to the sensor with wavelength 2 = 40 um. The thickness of the piezoelectric substrate (h1) is 
400 um, while the thickness of the magnetic sensitive layer (h3) varies to determine the optimal working 
point of the sensor. Subsequently, the simulation changes the thickness (h1) of the piezoelectric layer to 
determine the limit of the sensor’s working point when the wavelength à is constant based on the found 
optimal point (h3). 


Figure 1. The SAW-MO sensor structure 


2.2. Mathematical model for FEM simulation 

In this work, we set up a mathematical model of the SAW-MO sensor in a piezoelectric material 
environment, representing the forward and reverse piezoelectric effects with a set of electromechanical (1) 
and (2) [24], [25]. 


T = cS- eE (1) 
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D = eTS + Æ (2) 


i 
‘Lu 


Figure 2. The SAW-MO sensor structure parameters 


Table 1. Structural parameters of the sensor 


Parameters Symbol Value (um) 

Number of IDT finger pair n 25 pairs 
Width of each finger d Vary 
Left edge distance L2 2,000 
Right edge distance L1 2,000 
Length of IDT L— idt 990 
Length of magnetic sensitive layer Ls 4,960 
Thickness of piezoelectric layer h1 Vary 
Thickness of magnetic sensitive layer h3 Vary 
Gap between IDT and sensitive layer L3 20 
Sensor length L 10,980 


These equations involve parameters such as the elastic matrix (c), piezoelectric matrix (e), and 
dielectric constant matrix (€) along with stress (T), strain (S), electric displacement (D), and electric field (E) 
(see in Table 2). We also use motion (3) and (4) [24] to describe particle movement in the solid, with phase 
velocity, particle displacement, and potential as variables. 


ð2 
V-T =p (3) 
E = -V® (4) 


Where u is the particle displacement in the solid, Ø is the potential, and p is the density of the piezoelectric 
material. Here the operators are: 


0 0 0 
“lie 5a 3e (5) 
ð 0 
TA 0 0 0 TA TA 
ð 0 
pat o 2 So (6) 
0 0 0 
0 0 Ia TE T 0 


Combining (3) to (6) with (1), (2) and the condition V - D = 0, we can express the constitutive (7) in 
index form for the AIN piezoelectric material [4]. 


d7ut d7uf d?o 
p i_¢a k SS 
A at? ijkl axax; KU dxyax; 
duk a?o (7) 
jkl dxıdx jk dxęķdxj 


This equation involves the density, elastic coefficient, piezoelectric coefficient, mechanical 
displacement, and electric field of AIN, as well as indices i, j, k, l=1, 2, 3. Similarly, we have the constitutive 
(8) for FeNi magnetic sensitive material, which does not involve electric or displacement components [4]. 
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The wave used in our study is the Rayleigh wave (Vp), and the velocity relationship between this 
wave and the shear wave (V;) is described in (9) and (10) [26], [27], with anisotropy coefficient A and shear 
wave velocity V; calculated differently for isotropic and anisotropic materials [24], [26]. 


VR _ 0.436+c12/C11 


Vs 0.5+C42/C14 (9) 
VR [24 eee 1/2 
ve Fa | (10) 


We note that Vs is always greater than the surface acoustic wave velocity. 


Table 2. Physical parameter set of AIN, Al, and FeNi materials [11], [28]-[30] 
Value 

AIN Al FeNi 

Elastic stiffness coefficients (10''N.m?) Ch 4.10 


Material parameters 


Cr 1.49 
Ci3 0.99 
C33 3.89 
Cu 1.25 
C66 1.25 
Piezoelectric coefficients (C.m® €31 €33 0.58 
C24 1.55 
0.48 
Dielectric constants E 9.0 
€22 9.0 
€33 9.0 
Density (kg.m™) P 3,300 2,697 8,380 
Young’s modulus (GPa) E 70.3 E-H 
Poison ratio v 0.345 0.3 


Based on the sensor structure illustrated in Figure 2, the surface acoustic wave propagates through 
the delay-line element, which is coated with a magnetic sensitive layer made of FeNi with a variable 
thickness. Consequently, the surface acoustic wave experiences two factors. First, the mass (m) of the FeNi 
magnetic sensitive layer reduces the surface acoustic wave velocity [11], which can be calculated using (12). 
Second, the surface acoustic wave velocity gradually increases and approaches the shear velocity when the 
thickness of the FeNi sensitive layer increases [24], [26] or the AIN piezoelectric layer decreases. The mass 
is calculated as (11): 


m = p,x(h3xLsxDs) (11) 

In this simulation, the depth (Ds) and the length (Ls) of FeNi sensitive layer are constant, while its 
thickness (h3) varies, and its density (pp) is defined as the constant k = Ls X Ds. Therefore, the mass of the 
FeNi magnetic sensitive layer can be expressed using: 


m = p,xkxh3 (12) 


Moreover, the elastic coefficients of the FeNi magnetic sensitive material layer are determined by the 
relationship [11]. 


cE = E(1-vp) 
11 (14+ vp)(1-2 vp) 
es ee 
i2 = (1+vF)(1-2vF) (13) 
F _ E(Q-2vp) 
CAA m NA 


2(1+vp)(1-2 ve) 
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Where E is the Young’s modulus defined by the E-H characteristic and the poison coefficient (vp) is equal to 
0.3. H is the measured magnetic field intensity. 

To simulate the sensor’s working response accurately, the E-H characteristic described in [11] is 
used to determine the sensor’s working points. Specifically, as the measured magnetic field intensity (H) 
changes, it alters the young’s modulus (E) of the magnetic sensitive material layer, which affects the 
propagation velocity of surface acoustic wave (Vp). Combining this with the relation AVp/Vp = Af /f0 [31], 
where fo is the center frequency of the sensor, determines the frequency change of the sensor. By varying the 
thickness (h3) of the FeNi sensitive layer at H = 0 Oe, we can establish the relationship between the center 
frequency and the thickness of the FeNi sensitive layer to identify the optimal working point, which 
maximizes the acoustic wave velocity or resonant frequency at a given thickness. Subsequently, additional 
simulations are performed around the optimal working point to determine the sensor’s working response and 
draw conclusions regarding the best sensitivity of the magnetic sensor. With the optimal working point fixed, 
the thickness of the AIN piezoelectric layer is varied to determine the central frequency limit point of the 
sensor with a constant wavelength. This serves as the basis for increasing the center frequency of the sensor, 
which hopefully enhances the sensitivity of the sensor when working at lower measuring ranges. 


2.3. FEM simulation process 

To simulate the sensor’s working response, we used the constitutive (7) and (8) and data from Table 2. 
We performed FEM simulation in the following steps: firstly, entering the physical parameters (Table 2) and 
declaring the structure parameters (Table 1); secondly, building the sensor structure (as shown in Figure 2), 
meshing and changing material layers; thirdly, applying voltage to the IDT-in and setting boundary 
conditions; finally, solving and reading results. The simulation process is illustrated in Figure 3. 


Declare material, structure construction, meshing 
and material exchange. 


Apply the voltage on IDT-in and set boundary conditions. 


Simulation ends. 


Figure 3. FEM simulation process images 


We ensured accuracy while improving simulation speed by meshing the sensor with thicker meshing 
density at the top where the IDTs and magnetic sensitive layer FeNi are located and thinner density at the 
bottom. With a meshing step of about 2/40 = 1 um, the model had around 113,556 elements, ensuring 
convergent results. The boundary conditions were set as follows: the top surface had a free mechanical 
condition and zero voltage, the bottom had a fixed mechanical condition and grounding voltage, and the 
junction between two materials had free mechanical conditions and continuous voltage. All side edges of the 
sensor had free conditions. 

We applied a square pulse voltage of 100 V and pulse width of 10 ns to the IDT-in, then let the FEM 
calculation run long enough for the surface acoustic wave to propagate from IDT-in to IDT-out (about 
1,500 ns), with a signal sampling period of 0.3 ns to satisfy Shannon’s theorem. The output voltage was read 
on the IDT-out in the time domain and converted to the frequency domain using the FFT spectral density 
transform algorithm to determine the resonant frequency (f,) and surface acoustic wave velocity. We 
iteratively repeated this process to find the optimal working point by changing the FeNi sensitive layer 
thickness and to determine the working responses around the optimal point as the measured magnetic field 
intensity (H) changes. We also changed the AIN piezoelectric layer thickness to determine the limit point of 
the center frequency and the sensor’s working response when the center frequency is increased, with the aim 
of improving the sensor’s sensitivity when measuring at a lower range. 
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3. RESULTS AND DISCUSSION 

The simulation results presented in Table 3, Figures 4 and 5 reveal a two-stage behavior in the 
characteristic. In the first stage, increasing the thickness of the FeNi magnetic sensitive layer up to 1,060 nm 
results in an increase in the center frequency and surface acoustic wave velocity. However, in the second 
stage, further increases in the thickness of the FeNi layer lead to a decrease in the center frequency and 
surface acoustic wave velocity. These results are consistent with the theoretical framework proposed in this 
study and corroborated by previous research [4], [11], [15]. 


Table 3. Center frequency (fọ) and surface acoustic wave velocity (Vp) at the different thickness points of the 
FeNi magnetic sensitive layer (h3), Ve = Ax fo, h1 = 400 um and A = 40 um 
h3 (nm) fo (Hz) Vg (m/s) 

340 139,518,230 5,580.7 
420 139,587,400 5,583.5 
580 139,766,440 5,590.7 
660 139,933,270 5,597.3 
760 140,051,270 5,602.1 
820 140,157,060 5,606.3 
920 140,238,440 5,609.5 
1,000 140,303,550 5,612.1 
1,060 140,384,930 5,615.4 
1,140 140,319,820 5,612.8 
1,210 140,177,410 5,607.1 
1,280 139,672,850 5,586.9 
1,340 139,229,330 5,569.2 
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Figure 4. Relationship between the fọ and FeNi layer Figure 5. Relationship between the Vp and FeNi 
thickness (h3) at H = 0 Oe layer thickness (h3) at H = 0 Oe 


Based on the simulation results, we identify the optimal working point of the sensor at h3 = 1,060 nm, 
when h1 = 400 um and A = 40 um remain unchanged. This conclusion follows from the fact that the center 
frequency and surface acoustic wave velocity reach their maximum values at this thickness. To gain a clearer 
understanding of the process used to determine the center frequency of the sensor at different thicknesses 
(h3) of the FeNi magnetic sensitive layer, Figure 6 displays the frequency response of the sensor at three 
thickness points: 1,000 nm, 1,060 nm, and 1,210 nm. We determine the center frequency of the sensor for 
each thickness point by identifying the resonance peak of the frequency response. 

We simulated the sensor’s working response at different FeNi magnetic sensitive layer thicknesses, 
including 1,000 nm, 1,060 nm, and 1,210 nm, around the optimal working point. The simulation results are 
presented in Figure 7 and Tables 4 to 6. At h3=1,000 nm, the sensor’s response range is H = [0 + 109] Oe, 
with corresponding output frequencies of f = [140.30355 + 139.35954] MHz, or frequency shift 
Af = [0 + — 944.01] kHz. The sensor has a center frequency of fo = 140.30355 MHz when H = 0 Oe. 
Similarly, at h3 = 1,060 nm, the sensor’s response range is H = [0 + 89.0] Oe, with corresponding output 
frequencies of f = [140.38493 + 139.4694] MHz, or frequency shift Af = [0 + — 915.53] kHz. The sensor 
has a center frequency of fọ = 140.38493 MHz when H = 0 Oe. At h3 = 1,210 nm, the sensor’s response 
range is H = [0 + 112.0] Oe, with corresponding output frequencies of f = [140.17741 + 139.12269] MHz, 
or frequency shift Af = [0 + —105.472] kHz. The sensor has a center frequency of fo = 140.17741 MHz 
when H = 0 Oe. 
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Figure 7. Working response of the sensor at the 
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Table 4. Frequency response (f) and frequency shift (Af) to measured magnetic field intensity (H) of the 
sensor at h3 = 1,000 nm, h1 = 400 um and A = 40 um 


HO) f (Hz) AfD 
0.0 140,303,550 0 
15.0 140,295,410 -8,140 
30.7 140,250,650 -52,900 
40.0 140,205,890 -97,660 
48.0 140,140,790 -162,760 
61.0 139,937,340 -366,210 
81.0 139,542,640 -760,910 
89.0 139,432,780 -870,770 
100.0 139,367,680 -935,870 
109.0 139,359,540 -944,010 


Table 5. Frequency response (f) and frequency shift (Af) to measured magnetic field intensity (H) of the 
sensor at h3 = 1,060 nm, h1 = 400 um and A = 40 um 


H (Oe) f (Hz) Af (Hz) 
0.0 140,384,930 0 
15.0 140,368,650 -16,280 
30.7 140,291,340 -93,590 
40.0 140,193,680 -191,250 
48.0 140,071,610 -313,320 
61.0 139,798,990 -585,940 
81.0 139,522,300 -862,630 

89.0 139,469,400 -915,530 


Table 6. Frequency response (f) and frequency shift (Af) to measured magnetic field intensity (H) of the 
sensor at h3 = 1,210 nm, h1 = 400 um, and A = 40 um 


H (Oe) f (Hz) Af (Hz) 
0.0 140,177,410 0 
15.0 140,148,930 -28,480 
30.7 140,055,340 -122,070 
40.0 139,973,960 -203,450 
48.0 139,892,580 -284,830 
61.0 139,721,680 -455,730 
81.0 139,444,990 -732,420 
89.0 139,318,850 -858,560 

100.0 139,192,710 -984,700 
109.0 139,127,600 -1,049,810 
112.5 139,122,690 -1,054,720 


We determined the sensor’s sensitivity in two cases based on its working responses at different h3 
points. In the first case, the sensitivity is calculated over the full response scale of the sensor. At h3 = 1,000 nm, 
the response range H = [0 + 109] Oe, with an achieved sensitivity Sn of 8.661 kHz/Oe; at h3 = 1,060 nm, 
the response range H = [0 + 89] Oe, with an achieved sensitivity Sn of 10.287 kHz/Oe; and at h3 = 1,210 nm, 
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the response range H = [0 + 112.5] Oe, with an achieved sensitivity Sn of 9.375 kHz/Oe. The details are 
presented in Table 7 and Figure 8. In the second case, the sensitivity is calculated on the same response scale 
of the sensor H = [0 + 89] Oe. At h3 = 1,000 nm, the achieved sensitivity Sn is 8.661 kHz/Oe; at 
h3 = 1,060 nm, the achieved sensitivity Sn is 10.287 kHz/Oe; and at h3 = 1,210 nm, the achieved 
sensitivity Sn is 9.375 kHz/Oe. 


Table 7. The sensitivity (Sn) of sensor at points h3 when A = 40 um and h1 = 400 um 


The full scale < A 
h3 (nm) Sn (kHz/Oe) h3 (nm) Sn (kHz/Oe) 
0 0 0 0 
1,000 8.661 1,000 9.784 
1,060 10.287 1,060 10.287 
1,210 9.375 1,210 9.647 


In the second case, the sensitivity of the sensor is calculated on the same response scale of the 
sensor H = [0 + 89] (Oe), with achieved sensitivities of 8.661 kHz/Oe at h3 = 1,000 nm, 10.287 kHz/Oe at 
h3 = 1,060 nm, and 9.375 kHz/Oe at h3 = 1,210 nm. Details are shown in Table 7 and Figure 9. These 
results indicate that at h3 = 1,060 nm, the sensor exhibits the highest surface acoustic wave velocity and 
center frequency, resulting in the maximum sensitivity of 10.287 kHz/Oe in both cases. Our findings are in 
line with previous research in [6], which shows that higher surface acoustic wave velocity or center 
frequency leads to greater sensitivity. Additionally, to identify the points on the response range of the sensor 
when the magnetic field intensity changes from 0 to 89 Oe, we examined the resonance peaks on the 
frequency responses, as shown in Figures 10 to 12. 
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We conducted simulations to investigate the influence of piezoelectric substrate thickness on the 
central frequency limit of the sensor. The results indicate that decreasing the piezoelectric substrate thickness 
increases the center frequency, as shown in Figures 13 and 14. This finding is consistent with previous 
studies in [6], [24], [26], [27], which suggest that as the piezoelectric substrate thickness decreases, the 
surface acoustic wave velocity approaches the shear wave velocity. 
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Figure 12. Frequency response of FeNi/IDT/AIN Figure 13. Center frequency response with varying h1 at 
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Figure 14. Frequency response of the sensor at h1 with keeping h3 = 1,060 nm and 4 = 40 um 


However, when the piezoelectric substrate thickness is reduced to less than one wavelength (see in 
Figure 14), the frequency characteristic at h1 = 35 um no longer exhibits a resonance peak, leading to the 
indeterminability of the center frequency, unlike points with greater thickness. This phenomenon can be 
explained by the fact that when voltage is applied to the IDT-in, the piezoelectric material undergoes 
deformation due to the reverse piezoelectric effect from the surface to a depth of approximately one 
wavelength (A), as stated in [26]. Thus, if the piezoelectric substrate thickness is greater than approximately 
one wavelength, the deformation region on the piezoelectric substrate is preserved. However, if the thickness 
is less than one wavelength, the deformation region is not preserved, causing the surface acoustic wave to be 
cut off at the lower peak, and as a result, the frequency characteristics are not preserved. 

To preserve the frequency characteristics of the sensor, we can reduce the wavelength to A = 30 um 
or increase the center frequency. The simulation results, which are shown in Figure 15, demonstrate that the 
frequency center fọ = 218.46517 MHz when the wavelength is reduced to A = 30 um. The working range 
of the sensor is [0 + 89] Oe, as shown in Figure 10, and the highest sensitivity is achieved when the 
magnetic sensitive layer thickness h3 is 1,060 nm. However, the sensitivity is lower in the segment from 0 to 
35 Oe, as seen from the frequency shift characteristic of the sensor at h3 = 1,060 nm. To increase the 
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sensitivity of this segment, we can reduce the piezoelectric substrate thickness, which also reduces the 
fabrication and simulation time, and increase the center frequency of the sensor. The response simulation 
results of the sensor are shown in Figure 16. The simulation results show that the center frequency 
fo = 218.46517 MHz, the working range reaches from 0 to 33.1 Oe with the frequency shift of 248.21 kHz 
when the sensor has the parameters h1 = 35 um, A = 30 um, and h3 = 1,060 nm. On the other hand, the 
structure of the sensor with the parameters h1 = 400 um, 2 = 40 um, and h3 = 1,060 nm only gives the 
frequency shift of 118.79 kHz when considering the working range from 0 to 33.1 Oe. These results 
demonstrate that increasing the center frequency increases the sensitivity of the sensor, as shown in Table 8. 


I x I H 
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Figure 15. Frequency response of the sensor when h1 = 35 um, à =30 um, h3 = 1,060 nm 
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Figure 16. Frequency shift response of sensor at h3 = 1,060 nm 


Table 8. The Sn and Af of two structures with working range from 0 to 33.1 Oe 
Structure parameters 
him) Am) h3 (nm) Af (kHz) Sn (kHz/Oe) 
400 40 1,060 118.79 3.589 
35 30 1,060 248.21 7.499 
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Obviously, in this research, we have optimized the design parameters of the FeNi/IDT/AIN sensor 
structure with the sensitivity of 10.287 kHz/Oe which increases ~1.183 times compared to the sensitivity 
(8.698 kHz/Oe) of the FeNi/IDT/AIN sensor structure and ~4.387 times compared to the sensitivity 
(2.345 kHz/Oe) of the FeNi/IDT/LiNbO3 sensor structure in previous work [6]. This increasing sensitivity is 
significant and indicates that it is necessary to perform simulations to optimize the sensor parameters. 


4. CONCLUSION 

In this paper, we used finite element method simulations to investigate the SAW-MO sensor with a 
FeNi/IDT/AIN structure. Firstly, we explored the optimal thickness (h3) of the FeNi magnetic sensitive layer 
to achieve the best sensor performance. Next, we examined the impact of the magnetic sensitive layer 
thickness on the sensor’s sensitivity. We also considered the influence of the AIN piezoelectric substrate 
thickness on the sensor’s center frequency. Our results demonstrate that the optimal structure of the sensor is 
achieved at h3 = 1,060 nm, with a highest surface acoustic wave velocity of 5,615.4 m/s and a highest 
center frequency of 140.38493 MHz. We observed that the sensitivity of the sensor reaches its maximum 
value of 10.287 kHz/Oe at h3 = 1,060 nm. Additionally, we found that reducing the thickness of the AIN 
piezoelectric substrate to h1 = 35 um at the optimal point h3 = 1,060 nm resulted in an increased 
sensitivity of 7.499 kHz/Oe and a working range from 0 to 33.1 Oe. In comparison, the sensitivity of the 
sensor structure with h1 = 400 um, A = 40 um, and h3 = 1,060 nm was only 3.589 kHz/Oe within the 
same working range. In the future, the effect of the surface morphology of magnetic sensitive layer on the 
sensor’ sensitivity will be conducted and experiments will be performed. 
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